syms N11 N21 N31 N41 N12 N22 N32 N42 N13 N23 N33 N43 N14 N24 N34 N44 ...
     xD11 xD21 xD31 xD41 xD12 xD22 xD32 xD42 xD13 xD23 xD33 xD43 ...
    xD14 xD12 xD34 xD44 yD11 yD21 yD31 yD41 yD12 yD22 yD32 yD42 ...
    yD13 yD23 yD33 yD43 yD14 yD12 yD34 yD44

B = [N11 N21 N31 N41;
    N12 N22 N32 N42;
    N13 N23 N33 N43;
    N14 N24 N34 N44];

D = [xD11 xD21 xD31 xD41;
    xD12 xD22 xD32 xD42;
    xD13 xD23 xD33 xD43;
    xD14 xD12 xD34 xD44;
    yD11 yD21 yD31 yD41;
    yD12 yD22 yD32 yD42;
    yD13 yD23 yD33 yD43;
    yD14 yD12 yD34 yD44];